Discovery of >350 GeV Gamma Rays from the BL Lacertae Object lES 

2344+514 
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ABSTRACT 

We present the discovery of >350 GeV gamma-ray emission from the BL Lacertae 
(BL Lac) object lES 2344+514 with the Whipple Observatory 10m gamma-ray 
1^ , telescope. This is the third BL Lac object detected at very high energies (VHE, 

(SJ ! £; > 300 GeV), the other two being Markarian 421 (Mrk 421) and Mrk 501. These 

three active galactic nuclei are all X-ray selected and have the lowest known redshifts of 
any BL Lac objects currently identified. The evidence for emission from lES 2344+514 
■ comes mostly from an apparent flare on 1995 December 20 (universal date) during 

which a Qa excess was detected with an average flux of I(>350 GeV) = 6.6 ± 1.9 x 10^^^ 
Q^. photons cm~^ s^^. This is approximately 63% of the VHE emission from the Crab 

Q , Nebula, the standard candle in this field. Observations taken between 1995 October 
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and 1996 January, excluding the night of the flare, yield a 4cr detection indicating a 
flux level of I(>350 GeV) = 1.1 ± 0.4 x IQ-^^ photons cm'^ s-\ or about 11% of the 
VHE Crab Nebula flux. Observations taken between 1996 September and 1997 January 
on this object did not yield a signiflcant detection of a steady flux, nor any evidence 
of flaring activity. The 99.9% confidence level upper limit from these observations is 
I(>350 GeV) < 8.2 x lO^^^ photons cm^^ g-i, ^8% of the Crab Nebula fiux. The low 
baseline emission level and variation in nightly and yearly flux of lES 2344+514 are 
the same as the VHE emission characteristics of Mrk 421 and Mrk 501. 



Subject headings: BL Lacertae objects: individual (lES 2344+514) — gamma rays: 
observations 



1. Introduction 



The Whipple collaboration has discovered the two extragalactic sources of very high energy 
(VHE, E^300 GeV) gamma rays Markarian 421 ( [Punch et al. 1992D and Markarian 501 ( |Quinn 
|et al. 1996| ). Both of these are BL Lacertae (BL Lac) objects, a sub-class of the blazar class 
of active galactic nucleus (AGN). Blazars are the only class of AGN detected above 100 MeV 
by the Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma-Ray 
Observatory (CGRO) ([Fichtel et al. 1994| ; [Thompson et al. 1995| ; [Thompson et al. 1996| ) and BL 
Lac objects make up a signiflcant fraction of the blazars detected by EGRET. The most striking 
feature of the VHE emission from Markarian (Mrk) 421 and Mrk 501 is the day-scale or shorter 
flaring which has produced emission as high as 10 times the Crab flux ( Gaidos et al. 1996| ), but 



flares on the order of 1/2 the Crab flux are more typical (Buckley et al. 1996; Quinn et al. 1996; 



Quinn et al. 1997; McEnery et al. 1997). The baseline emission level for these two objects can be 
very low; Mrk 501 was initially detected to have a VHE flux 8% that of the Crab Nebula ( Puinn| 



et al. 1996[ ) and Mrk 421 has dropped below the detection limit of the Whipple Observatory 
gamma-ray telescope for as long month ( [Buckley et al. 1996[ ; [McEnery et al. 1997[ ). 



Initial searches for extragalactic sources of VHE gamma rays by the Whipple collaboration 
concentrated mainly on blazars detected by EGRET or AGN of any type with small redshift 



( Kerrick et al. 1995 ). This approach led to the detection of Mrk 421, but yielded only upper 
limits for the other objects studied. Beginning in the spring of 1995, the Whipple collaboration 
initiated a more focussed observation campaign on nearby BL Lac objects. The search list was 
mainly drawn from the work of Perlman et al. (1996), who identifled BL Lac objects from the 
Einstein Slew Survey sources ( Elvis et al. 199^ ) f] This list contained all of the prominent BL Lac 
objects (e.g., BL Lacertae, PKS 2155-304, Mrk 421), but also signiflcantly increased the number 



^We note that Mrk 501 was chosen as a candidate independently because of its similarities to Mrk 421. 



- 3 - 



of BL Lac objects with known redshifts. In addition, because these objects were drawn from an 
X-ray survey, they were more hkely to have emission spectra similar to Mrk 421 and Mrk 501. We 
hmited the search to those sources with low redshifts (initially z < 0.1, but eventually extended to 
z < 0.2) in order to reduce the chances that the VHE emission would be significantly attenuated 
by interaction with the background intergalactic infrared radiation fields ( pould &: Schreder 



1967 ; Stecker, de Jager, fc Salamon 1992| ). The first success of this observing program was the 



detection of Mrk 501 as a gamma-ray source. As Mrk 501 was not identified in EGRET catalogs 
as a significant source, it highlighted the ability of ground-based gamma-ray telescopes to not 
only complement the results of the space-based gamma-ray telescopes, but to augment them. We 
report here on the second source detected in this BL Lac object search program, lES 2344+514. 



lES 2344+514 was only recently identified as a BL Lac object (Perlman et al. 1996), based 
on its lack of optical emission lines with observed equivalent width greater than 5A and its Ca 
II "break strength" being smaller than 25%. The latter criteria is indicative of the presence of a 
power law continuum and the former eliminates quasars. Perlman et al. (1996) determined this 
object to have a redshift of z = 0.044 based on absorption lines; it had no evident emission lines. 
This makes it the third closest known BL Lac object, after Mrk 421 and Mrk 501. Perlman et al. 
(1996) derive a 2 keV X-ray flux of 1.14 /uJy, roughly 1/3 the flux detected for Mrk 421 and Mrk 
501 in that same work, and measured an optical magnitude of my = 15.5 with no galaxy light 
subtraction. Measurements taken with the Very Large Array radio interferometer indicate that 
its radio emission is "point-like," with more than 80% of its flux being from an unresolved point 
source (patnaik et al. 1992| ; [Perlman et al. 1996| ). The Gr een Bank radio survey lists its 5 GHz 



radio flux as 231+25 mJy, which is about 1/3 and 1/4 the 5 GHz flux of Mrk 421 and Mrk 501, 
respectively. 

In the following sections we outline the observation and analysis techniques (§ ^), detail the 
observations taken by the Whipple collaboration on lES 2344+514 between 1995 October and 
1997 January (§ |3|), and flnally discuss the outcome and implications of those observations (§ 



2. Observation and Analysis Techniques 

The VHE observations reported in this paper were made with the atmospheric Cherenkov 
imaging technique ( Pawley fc Weekes 1995 ; Reynolds et al. 1993] ) using the 10 m optical reflector 



located at the Whipple Observatory on Mt. Hopkins in Arizona (elevation 2.3 km) (Cawley et 



al. 19901) . A camera, consisting of photomultiplier tubes (PMTs) mounted in the focal plane of 



the reflector, records images of atmospheric Cherenkov radiation from air showers produced by 
gamma rays and cosmic rays. For most of the observations reported here, the camera consisted of 
109 PMTs (each viewing a circular fleld of 0?259 radius) with a total fleld of view (FOV) of 3°. In 
1996 December, 42 additional PMTs were added to the camera, increasing the field of view to 3?3. 
Because only a small fraction of the data presented here was taken with more than 109 PMTs, the 
analysis of this data only uses the original 109 PMTs. This makes the results consistent with the 
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rest of the 1996/97 season and allows one set of parameter cuts to be used in the analysis. The 
telescope is sensitive to gamma rays with energies between approximately 200 GeV and 20 TeV. 

We characterize the Cherenkov images' roughly elliptical shapes and locations and orientations 
within the telescope FOV using a moment analysis (Reynolds et al. 1993). The gamma-ray 



selection applied here utilizes the Supercuts95 criteria (Table 1; cf. Quinn et al. 1996; Catanese et 



al. 1996), which were optimized on contemporaneous Crab Nebula data to give the best sensitivity 
to point sources. Monte Carlo simulations indicate that this analysis results in an energy threshold 
of -^350 GeV and an effective area of 3.5 x 10^ cm^. Details of the methods used to estimate the 



energy threshold and effective area are given elsewhere ([Mohanty et al. 1997| . The large effective 
area makes this technique very sensitive to short-term variability in sources. 



2.1. On/OfF Observations 

The traditional mode of observing potential sources with the Whipple Observatory gamma-ray 
telescope, and the method which is still usually used to confirm source signals and derive source 
spectra is the On/Off mode of observation. In this type of observation, the putative source is 
tracked continuously for 28 minutes with the center of the telescope FOV positioned at the source 
location. In a 28 minute control, or "Off," run the telescope tracks a position offset by 30 minutes 
in right ascension but with the same declination as the putative source. This Off-source run begins 
exactly 30 sidereal minutes before or after the start of the On-source run so that the telescope 
tracks the same range of elevation and azimuth for both runs. Night sky brightness differences 
between the two observing fields are equalized off-line by software padding ( Cawley et al. 19"93|) . 



The significance of any excesses or deficits in the observations are calculated using the maximum 
likelihood method of Li & Ma (1983). 

This data collection technique has been shown to produce statistically stable gamma-ray 
count rates from the Crab Nebula over the course of several months ( iQuinn 1997 ), and to allow 



consistent, reproducible spectra to be derived for the existing sources ( Mohanty et al. 1997 ; Hillas 



et al. 1997| ; Zweerink et al. 1997). However, this method requires excellent weather conditions to 



obtain consistent results between the On-source and Off-source runs so it reduces the observing 
time available. Also, because it only uses a small fraction of its gamma-ray like events to estimate 
the background, it does not give the most statistically accurate background estimates possible. 



2.2. Tracking Observations 

For the reasons listed above, the more frequent mode of observation, particularly when 
searching for new sources, is the Tracking mode. This mode of data collection for the Whipple 
Observatory gamma-ray telescope has been described previously ( Kerrick et al. 1995| ; Puinn e"^ 
|al. 1996 ), but we do so again here for clarity and because some differences in the method of 
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significance calculation are used. In this mode, only the On-source position is tracked, in runs of 
28 minute duration, and no control observations are taken. To estimate the expected background, 
we use those events which pass all of the Supercuts95 criteria except orientation (characterized by 
the a parameter). We use events with values of a between 20° and 65° as the background region, 
as indicated in Figure |l|. The lower bound gives a 5° buffer between the On-source region and the 
Off-source region, to allow some room for spillover events from the source region. The upper bound 
is chosen because the region beyond 65° is less stable from run to run and source to source due to 
edge of field effects. As shown in Appendix this non-standard method of background estimation 
gives reliable results for a variety of observation fields. Using contemporaneous non-source data, 
we find that the factor (called a tracking ratio) which converts the number of events in the region 
a = 20° — 65° to an estimate of the On-source background (q = 0° — 15°) is r = 0.292±0.004 for 
the 1995/96 season and r = 0.316±0.004 for the 1996/97 season. 

To convert On-source and background counts to a significance, S, in the tracking analysis, we 
use a simple propagation of errors formula: 

S= , ^on-r*N,,, 

^Non + r2 * Nhkd + iV^d * 

where Non is the number of events in the source region (a < 15°), Nbkd is the number of events 
in the background region (a = 20° — 65°), and r it Ar is the tracking ratio and its statistical 
uncertainty. The maximum likelihood method of Li & Ma (1983) cannot be used when the ratio 
for converting the background counts to a background estimate has any uncertainty. Li &: Ma 
(1983) show that the estimate of the significance from equation (||) tends to be conservative when 
r < 1. 



2.3. Flux and flux upper limit estimation 

After Supercuts95 are applied to an On/Off or Tracking data set, we obtain a statistical 



significance as calculated with the methods described in § 2A and § |2.2| , respectively, and a 
corresponding gamma-ray count rate in terms of counts per minute. If the excess is not statistically 
significant, we calculate a 99.9% confidence level (CL) upper limit on the count rate using the 
method of Helene (1983). To convert these count rates or upper limits to fluxes, we first express 
them as a fraction of the Crab Nebula count rate for the same season. This corrects for season to 
season variations which affect the telescope response, and therefore its gamma-ray count rate. For 
the 1995/96 season, the gamma-ray count rate for the Crab Nebula with a Supercuts95 analysis 
is 1.58±0.05/min. In 1996/97, the Crab Nebula rate was 1.69±0.07/min. Analysis of the Crab 
Nebula data shows that for runs taken under good weather conditions, the gamma-ray count rate 



does not change significantly within a season but can change significantly between seasons ( Quinn 



1997. So, we can assume the gamma-ray count rate for a source can be reliably compared with 



that of the Crab Nebula over the course of a season. 
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Once we have the flux or flux limit expressed as a fraction of the Crab Nebula count rate, we 
multiply it by the Crab Nebula flux (in units of photons cm~^ s~^) above the energy threshold of 
the observations. We determined the energy threshold for the observations reported here to be 
350 GeV. The integral Crab Nebula photon flux is /(> 350 GeV) = (1.05 ± 0.24) x 10"^° photons 
cm~^ s^^ (Hillas et al. 1997). For flux estimates, we propagate the uncertainty in the Crab Nebula 



gamma-ray count rate and flux through to the final flux estimate, so the flux uncertainties have 
a significant contribution from the Crab Nebula flux uncertainty. The significances of an excess 
should thus not be estimated by the photon flux uncertainties. For the flux upper limits, we do 
not propagate through the Crab Nebula uncertainties. Upper limits are an estimate of the flux 
that could be present in the data set but not produce a significant excess. This is most accurately 
derived from the count rate because that is what determines the statistical significance of the 
excess. The Crab Nebula count rate and flux uncertainties affect only the normalization, so the 
flux upper limits quoted in terms of photons cm~^ s~^ have an uncertainty of ~25%, mainly from 
the uncertainty in the Crab Nebula photon flux. 

This flux calculation takes advantage of the fact that the VHE Crab Nebula flux is steady 
over at least a 7 year period (Hillas et al. 1997) so that changes in the Crab Nebula count rate 



are most likely due to changes in the telescope sensitivity or threshold. The one disadvantage of 
this method is that it assumes the spectrum of the putative source is similar to that of the Crab 
Nebula (oc E~'^'^). Preliminary estimates of the spectrum for lES 2344+514 (to be presented in 
future work) indicate that it is not significantly different than this. 

As shown below, lES 2344+514 is observed with an average elevation of ~ 55°, much lower 
than the average elevation of the Crab Nebula observations (~ 72°). While the telescope trigger 
rate does vary significantly with elevation, Supercuts95 analysis of the Crab Nebula indicates no 
significant variation in the gamma-ray count rate with elevation down to 50°. So, normalizing to 
the Crab Nebula count rate should be valid over that elevation range. 



3. Observations and Results 

For the 1995/96 observing season, lES 2344+514 was observed between October and January. 
After filtering the runs for bad weather and instrumental problems, the data set consists of 
19 Tracking runs (exposure = 9.0 hours) and 26 On/Off pairs (exposure = 11.5 hours). The 
elevation of these runs varied from 37° to 70° with a mean of 55°. For the On/Off pairs, using the 
Supercuts95 gamma-ray selection criteria, we obtain a significance of 2.9a with a corresponding 
excess count rate of 0.19±0.07/min or a 99.9% CL upper limit of <0.40/min. This converts to a 
flux upper limit of <0.25 times the VHE Crab Nebula flux rate or an integral flux upper limit of 
I(>350 GeV) < 2.6 x IQ-^^ photons cm'^ s'^ 

Combining all of the On-source runs whether taken as On/Off pairs or Tracking runs, and 



using the Tracking analysis described in § 2.2, we obtain a significance of 5.8a and count rate of 
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0.25ib0.04/miii. This translates into a flux of 0.16ib0.03 times the Crab Nebula flux or I(> 350 
GeV) = 1.7 lb 0.5 X 10^^^ photons cm~^ s~^. The higher significance is due to a combination of 
having more data and the tracking analysis having a more statistically accurate estimate of the 
background. The count rates for the two data sets are consistent within statistical errors. These 
results are summarized in Table |2|. Figure Q shows the a distribution for the summed On-source 
runs, along with the a distribution for the non-source runs used to estimate the tracking ratio. 
The latter distribution is normalized to the predicted background level from the Tracking analysis. 

For the 1996/97 observing season, lES 2344+514 was observed between September and 
January. After filtering for bad weather and instrumental problems, the data consist of 38 On/Off 
pairs (exposure = 17.4 hours) and 16 Tracking runs (exposure = 7.5 hours). The runs were taken 
at elevations ranging from 46° to 70° with an average elevation of 64°. For the On/Off pairs, 
Supercuts95 analysis resulted in a -0.4(7 deficit, with corresponding 99.9% CL upper limit <0.15 
counts/min. This is equivalent to an upper limit on the flux of <0.09 times the Crab Nebula flux 
or I(>350 GeV) < 0.93 x 10^^^ photons cm^^ s^^. For the combination of the On-source runs 
taken as part of On/Off pairs or as Tracking runs, Supercuts95 analysis gives a 0.4cj excess with 
corresponding 99.9% CL upper limit <0.13 counts/min. This upper limit is equivalent to <0.08 
times the VHE Crab Nebula flux or I(>350 GeV) < 0.82 x 10"^^ photons cm'^ g-^. These results 
are summarized in Table ^. The a parameter distribution for the Tracking analysis is shown in 
Figure |2|. 



3.1. Variability 

As discussed in § |l], both Mrk 421 and Mrk 501 exhibit low baseline levels of VHE emission, 
but with comparatively high flux amplitude flares, typically lasting on the order of 1 day. So, this 
sort of behavior should be expected in the VHE emission from other BL Lac objects if they emit 
VHE gamma rays. Therefore, we have searched for flaring activity from lES 2344+514 in both the 
1995/96 and 1996/97 observing seasons. The light curves showing the daily excess or deflcit for all 
On-source data taken during the two seasons are shown in Figure ^. For 1996/97, the observations 
show no evidence of significant variability. The maximum daily excess has a significance of 1.8cr 
and the probability that the excesses and deficits are consistent with statistical fiuctuations 
about a common mean is 0.80. In contrast, for the 1995/96 season the probability that the 
observations are consistent with constant emission is 3.5 x 10~^. For comparison, the same analysis 
of the Crab Nebula data gives a x^ probability of constant emission of 0.78 in 1995/96 and 0.94 in 
1996/97. 

The most significant contribution to the variability in 1995/96 is 1995 December 20 (indicated 
by the arrow in Figure ^). On that night, 3 On/Off pairs and 1 Tracking run were taken on 
lES 2344+514. The On/Off pairs indicate a 4.2o- excess and a flux of 0.55+0.13 times the Crab 
Nebula flux, or I(>350 GeV) = (5.8 + 1.9) x 10^^^ photons cm~^ s~^. A Tracking analysis of all 
4 On-source runs reveals a 6. Oct excess with a corresponding flux of 0.63+0.11 times the Crab 
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Nebula flux, or I(>350 GeV) = (6.6 it 1.9) x 10~^^ photons cm~^ s~^. The increase in significance 
for the Tracking analysis is due to the additional run used. The fluxes for the On/OfF and Tracking 
analyses are consistent within errors. The a parameter distribution for the night of the flare is 
shown in Figure ^ clearly indicating the excess. There is no statistically significant variation in 
the flux within the night (see Figure |^ . 

The flare on 1995 December 20 constitutes slightly over 1/3 of the excess gamma rays detected 
from lES 2344+514 in the 1995/96 season. When the data from that night are removed, the 
Tracking analysis of all the remaining On-source runs reveals an excess of 4. Oct and corresponding 
flux of 0.11±0.3 times the Crab Nebula flux, or I(>350 GeV) = (1.1 ± 0.4) x 10"^^ photons cm-^ 
s~^. Thus, without the flare night, the evidence for emission from lES 2344+514 is marginal. 
The probability that the emission is constant about the mean for the non-flare nights is 0.003, 
indicating evidence for further variability at the 3a level. ^From Figure |3| it is evident that any 
additional flaring was not significant on day scales. Given the low significance of the overall excess, 
further study of this possible variability is not feasible. 



4. Discussion 

We have presented evidence of VHE emission from the BL Lac object lES 2344+514, as 
detected with the Whipple Observatory gamma-ray telescope. The emission is detectable only in 
1995/96 and only with unambiguous statistical significance on 1995 December 20. On the latter 
date, the flux was ~60% that of the Crab Nebula. The baseline emission appeared to be ~10% 
of the Crab Nebula flux during 1995/96 and was <8% in 1996/97. This low emission level with a 
high flux flare is exactly what would be expected from a VHE-emitting BL Lac object, based on 
what has been observed of the VHE emission from the other known extragalactic VHE sources, 
Mrk 421 and Mrk 501. In fact, the initial detection of Mrk 501 as reported by Quinn et al. (1996) 
is almost completely analogous to the observations reported here for lES 2344+514. The initial 
detection of Mrk 501 indicated a mean flux of 8% that of the Crab Nebula and in more than 
66 hours of observations, only a single flare was detected, with a flux equivalent to 50% of the 
Crab Nebula flux. The one difference is that in subsequent seasons, Mrk 501 has become brighter 
(Quinn et al. 1997). If it had become dimmer, it would not have been detected in subsequent 



seasons, just as happened with lES 2344+514 in 1996/97. 

The detection of lES 2344+514 is also consistent with Mrk 421 and Mrk 501 in that ah 
three are nearby. X-ray selected BL Lac objects (XBLs). By XBLs we mean here those BL Lac 
objects which have synchrotron spectra that extend into the UV to X-ray energy range. This 
extension makes them bright soft X-ray sources, hence the XBL label. The low redshift means its 
VHE emission will not be attenuated severely in the ~300 GeV energy range by interaction with 
background IR photons ( ^tecker, de Jager, fc Salamon 1996 ; Biller et al. 1997 ). BL Lac objects 



have been suggested as better candidates for VHE emission than the flat spectrum radio quasars 
(FSRQs) typically detected by EGRET ( Dermer &: Schlickeiser 1994| ) because BL Lac objects 
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lack broad emission lines in their optical spectra. The absence of such lines may indicate less 
VHE gamma-ray absorbing material at the source. XBLs should be the best blazar candidates 
for VHE emission if the gamma rays are produced from inverse Compton scattering of low energy 
photons by the same electrons which produce the synchrotron emission that dominates blazar 



spectra from radio to optical or X-ray wavelengths (e.g., Konigl 1981; Maraschi et al. 1992; Bloom 



fc Marscher 1996| ; ^ikora, Begelman, fc Rees 1994| ; ^ ermer, Schlickeiser, &: Mastichiadis 1992]). 



This is because the synchrotron spectrum of XBLs extends to higher energies than other blazars, 
like radio-selected BL Lac objects (RBLs) and FSRQs which have synchrotron cutoffs at optical 
wavelengths. This implies the presence of higher energy electrons in the jets of XBLs, so given 
similar magnetic fields and Doppler factors, XBLs should have an inverse Compton spectrum 
which extends to higher energies than those of RBLs or FSRQs (Sikora et al. 1994). This picture 
is consistent with the lack of VHE emission from nearby RBLs like BL Lacertae ( patanese et al. 



1997a) and W Comae ( Catanese et al. 1997b| ) which are EGRET sources (Catanese et al. 1997a; 



Thompson et al. 19961). However, it does not rule out models which produce gamma-ray emission 



as the product of extremely high energy protons in the AGN jet (e.g., Mannheim 1993). More 
detections of BL Lac objects will be needed to further investigate that issue. 

The spectral energy distribution (SED) of lES 2344+514 is indicated in Figure |^. Because the 
data are not contemporaneous, deriving limits on the jet parameters or making strong conclusions 
about the shape of the SED are not feasible. However, some general conclusions can be drawn. 
lES 2344+514 is dimmer at every wavelength than Mrk 421 and Mrk 501 and so it is perhaps not 
surprising that the VHE emission of lES 2344+514 appears to be weaker on average. If the X-ray 
power output were similar to that at 350 GeV, as it is in Mrk 421 (IMacomb et al. 1995| ; [Buckley et| 



al. 19961 ), the gamma-ray flux detected in the flare would require an X-ray flux at 2 keV of 12//Jy, 



an increase of a factor of 10 over its detection flux. This is an extreme, but not unprecedented, 
increase in X-ray flux. If the X-ray flux is higher than the gamma-ray flux, as it is with Mrk 
501 ( [Catanese et al. 1997c| ), the X-ray flux requirements are even greater. Contemporaneous 
observations of lES 2344+514 will do much to clarify the shape of its SED for comparison with 
the other VHE-emitting BL Lac objects. 

No other northern hemisphere imaging Cherenkov telescope has reported a detection or upper 
limit on lES 2344+514 at this time, though none were observing this object during 1995. Just 



as with Mrk 501, lES 2344+514 is not a source in the EGRET catalogs ([Fichtel et al. 1994[ ; 



Thompson et al. 1995; Thompson et al. 1996| ). The All-Sky Monitor on the Rossi X-ray Timing 



Experiment will allow better long-term multiwavelength monitoring of this object and continued 
monitoring of this object with Cherenkov gamma-ray telescopes will hopefully confirm this object 
as a VHE gamma-ray source in the near future. 
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A. Tracking mode observations 

In the Tracking observation mode the On-source and background data are recorded at identical 
times, so changes in external weather conditions or system-wide variations in the electronics 
should affect both the On-source and background regions similarly. The Tracking mode also allows 
more On-source data to be taken in a given amount of time because no separate Off-source run 
is needed. Finally, because more events (~3 times as many) are used in the background estimate 
than in the On/Off analysis, it can lead to more statistically accurate estimates of the mean 
background level. However, the Tracking analysis requires the distribution of background events 
versus a, the orientation parameter, to be the same from run to run in order to make accurate 
background estimates. Variations in the noise which do not affect the camera uniformly could 
change the a distribution. An example of such a variable noise source is the background light 
distribution due to stars within the FOV which changes with observation field and sky clarity. 
Also, a PMT which is turned off for a period of time, due to a bright star or malfunction, affects 
the telescope response in a non-uniform way. If the a distribution changes, the factor which 
converts the background events (i.e., those with a = 20° — 65°) to an On-source (a = 0° — 15°) 
background estimate changes as well. If this factor, called a tracking ratio, varies for different 
observations, the tracking analysis is not a viable means of estimating the background. 

To investigate these concerns, we have analyzed 148 hours of data taken on 20 non-source 
observation fields in the 1995/96 season and 99 hours of data taken on 19 non-source observation 
fields in the 1996/97 season. All runs were taken under good weather conditions. We analyzed 

these two data sets separately because they lead to different tracking ratios for the two seasons 
as shown below. These observations consist of Off-source runs and On-source runs for objects 
which do not reveal even a marginally significant signal. The only other selection made on these 
observation fields was that they not have any stars within the field of view which produce a 
current greater than ~100 /xA in more than one PMT. Stars which are that bright, such as ( 
Tauri 1?1 from the Crab Nebula, have a clear effect on the tracking ratio. The non-source fields 
used here have a range of mean background brightness levels due to low magnitude stars, and 
the elevations at which the observations were taken range from 20° to 85°. So, no particularly 
restrictive constraints are placed on the data that we use to determine a tracking ratio. This is 
as it should be; if we must be very restrictive about the type of field we use in a tracking ratio 
estimate, then the tracking analysis loses its usefulness. 

A tracking ratio of 0.292±0.004 is found for the 1995/96 data set and the 1996/97 data set 
indicates a tracking ratio of 0.316+0.004. The differences in the two tracking ratios are most 
likely due to the fact that the telescope camera PMT gains were set to different values in the two 
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seasons, resulting in a different response of the telescope to the Cherenkov and background light. 

To determine whether the different observation fields show evidence of inter-run, inter-source, 
or intra-season variability in the tracking ratio, we applied several tests. To investigate the 

stability of the tracking ratio from run-to-run, we calculated the probability that the tracking 
ratios for the individual runs were consistent with a single tracking ratio. For the 1995/96 data, 
the probability that the individual runs are drawn from the statistical distribution of a single 
tracking ratio is 0.23 and for the 1996/97 data, it is 0.16. Both are thus consistent with arising 
from a single tracking ratio. 

To test whether the different observation fields have different tracking ratios, we determined 
the significances of the excesses or deficits of the observations taken on the different fields and also 
calculated the probability that the objects' mean tracking ratios were drawn from a distribution 
with a common mean. If the tracking ratio varies with observation field, using a single tracking 
ratio will result in a significant excess or deficit in some fields. Among the 39 fields used in the 
tracking ratios for the two seasons, only one had a statistical excess or deficit of more than 2a. 
This is expected 86% of the time given the number of trials in this test. The probability that the 
means of the tracking ratios are drawn from a single mean is 0.73 for the 1995/96 data and 0.44 
for the 1996/97 season. So, both tests indicate that the tracking ratios for the individual fields are 
consistent with a single tracking ratio. 

The final test we applied was to see if the tracking ratio changed over the course of the 
season. To do this, we calculated the probability that the excess or deficit coTints per minute for 
the summed objects were consistent with statistical fluctuations about a constant mean. If the 
tracking ratio changes with time, it should produce a consistent excess or deficit that gives a poor 
fit to the constant mean hypothesis. The probability that the rate is constant is 0.68 for the 
1995/96 data set and 0.25 for the 1996/97 season. Thus, we have no evidence that the tracking 
ratio changes significantly over the course of an observing season. 

^Prom these tests, we conclude that the single tracking ratio can be used for objects which 
do not have particularly bright stars within the telescope's FOV. We do, however, see a significant 
difference in the tracking ratio between seasons. 
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Fig. 1. — The distribution of a parameter values for events passing all but the a parameter cut of 
Supercuts95. These data were taken on the Crab Nebula between 1996 October and 1997 January. 
The source and background regions used in the Tracking analysis described in the text are indicated. 



Table 1. Supercuts95 



Gamma-ray range 



maxi > 100 d.c.^ 
max2 > 80 d.c. 
size > 400 d.c. 
0?16 < length < 0?30 
0?073 < width < 0?15 
0?51 < distance < 1?10 
a < 15° 



^d.c. = digital counts. 
1.0 d.c. ft! 1.0 

photoelectrons. 
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Fig. 2. — Distributions of the a parameter for lES 2344+514 for events which pass all other 
Supercuts95 cuts. Shown are all On-source data taken (a) in 1995/96 and (b) in 1996/97. The 

dashed curves are the a distributions for the non-source runs used in each season to estimate a 
tracking ratio (see text for details) normalized to the predicted background in the a = 0° — 15° 
region. 
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Fig. 3. — Light curves for all On-source data taken in (a) 1995/96 and (b) 1996/97, for lES 

2344+514 passing Supercuts95. The uncertainties arc \a statistical errors. The arrow in (a) 
indicates the flare on 1995 December 20 (MJD 50071). In (b), MJD 50400 is 1996 November 13. 
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Fig. 4. — Distribution of the a parameter for all On-source data taken on lES 2344+514 on 1995 
December 20 for events which pass all other Supercuts95 cuts. On/Off data (a), and all On- 
source data (b) are shown. The dashed curve is the a distribution for the non-source runs used to 
estimate the 1995/96 tracking ratio (see text for details) normalized to the predicted background 
in the a = 0° — 15° region. 
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Fig. 5. — The hght curve for all On-source observations of lES 2344+514 taken on 1995 December 
20 for events passing Supercuts95. The error bars are la statistical uncertainties. 



Table 2. lES 2344+514 Observation Summary 







Exposure 


S 


Rate'' 


Flux*^ 


Flux'^ 


Epoch 


Data type 


(hours) 


(^) 


cts /min 


(Crab) 


(xlO"" cm~2 g-i^ 


1995/96 


On/Off 


11.5 


2.9 


<0.40 


<0.25 


<2.6 




Tracking'^ 


20.5 


5.8 


0.25+0.04 


0.16+0.03 


1.7+0.5 


1995 Dec 20 


On/Off 


1.38 


4.2 


0.87+0.21 


0.55+0.13 


5.8+1.9 




Tracking 


1.85 


6.0 


1.00+0.17 


0.63+0.11 


6.6+1.9 


1996/97 


On/Off 


17.4 


-0.4 


<0.15 


<0.09 


<0.93 




Tracking 


24.9 


0.4 


<0.13 


<0.08 


<0.82 


^The excess 


counts per 


minute, or 


the 99.C 


1% CL upper limit on 


that quantity, passing 



Supercuts95. 

^The flux, or 99.9% CL upper limit, is expressed as a fraction of the VHE Crab Nebula flux. 
^The integral flux, or 99.9% CL upper limit, is quoted for E > 350 GeV. 

'^In this table. Tracking refers to all On-source runs, including those taken as part of an 
On/Off pair. 
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Fig. 6. — The spectral energy distribution of lES 2344+514. The VHE observations from 1995/96 
(open square), 1996/97 (filled circle), and the flare flux of 1995 December 20 (filled star) are shown. 



Non-contemporaneous X-ray, optical and radio fluxes are indicated by the open circles (Perlman 
et al. 199^ ; Patnaik et al. 1992| ; Gregory fc Condon 1991] ). The optical flux is shown as an upper 
limit because the galaxy light contribution is not subtracted. In addition, the EGRET upper limit 
is shown (D.J. Thompson 1996, for the EGRET team, private communication). 



